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Abstract 

' In a previous paper, referred to as a "Mini-Landscape" search, we explored a 

■ "fertile patch" of the heterotic landscape based on a Kg-II orbifold with SO(IO) and 

Eg local GUT structures. In the present paper we extend this analysis. We find many 
models with the minimal supersymmetric standard model spectra and an exact R 
parity. In all of these models, the vector-like exotics decouple along D flat directions. 
We present two "benchmark" models which satisfy many of the constraints of a 
realistic supersymmetric model, including non-trivial Yukawa matrices for 3 families 
of quarks and leptons and Majorana neutrino masses for right-handed neutrinos with 
non-trivial See-Saw masses for the 3 light neutrinos. In an appendix we comment on 
the important issue of string selection rules and in particular the so-called "gamma- 
rule". 



1 Introduction 



The string landscape denotes the space of string vacua [1-3] . It is beheved that there are 
on the order of 10^°° or more possible vacuum configurations. However, at the moment, 
only a paltry few such vacua have properties looking anything like our own, i.e. with 
3 large space dimensions, the standard model gauge interactions and matter degrees of 
freedom, and a vanishingly small (in string units) cosmological constant. Much effort 
has gone into exploring the landscape in search of regions satisfying the latter requisite 
feature, while only a few groups have attempted to find the rest. If the string is to make 
contact with experiment, this situation must be inverted. 

Can string theory with 10^'^*' vacua make any predictions relevant for the LHC? 
It has been suggested that by exploring the entire string landscape one might obtain 
statistical data which could lead to probabilistic experimental statements [4,5]. Yet the 
clearest statement to-date is that standard-like models are exceedingly rare. This has 
been demonstrated in the context of orientifolds of Gepner models [6,7] and in the context 
of intersecting D-branes in an orientifold background [8-11]. Nevertheless they may be 
more prevalent in the heterotic string because of the simple fact that GUT groups, such 
as SO(IO) containing spinor representations, appear naturally. For explorations of the 
string landscape in the context of the free fermionic construction of the heterotic string, 
see [12-17]. However, even within the context of the heterotic string, standard model- 
like structure is not guaranteed. For example, it appears to be very difficult [18, 19], 
if not impossible [20,21], to find standard model-like structure in the heterotic string 
compactified on a Z3 orbifold. To summarize, standard-like models are very rare; quite 
possibly a negligibly small set in the entire landscape. If so, what might we conclude 
from this statistic? 

We suggest the following alternate strategy for obtaining low energy predictions from 
string theory. One should introduce some priors into ones statistical analysis. The first 
prior is that the theory has only 3 large space dimensions. The second is that the string 
vacuum includes the standard model. Perhaps within this subset one may find some 
statistical correlations which can be useful. Indeed, it is also important to verify that 
the standard model actually sits in the string landscape. Of course, in order to make 
this analysis tractable, one may need to include additional priors. In particular, one may 
require that the theory is supersymmetric at the string scale and that below the string 
scale the spectrum is that of the MSSM. Such theories typically have of order 100 moduli 
(geometric or others). The Yukawa and gauge couplings of the theory will generically be 
functions of these moduli. In the supersymmetric limit of the theory, one would hope to 
be able to tune the moduli in order to obtain acceptable low energy physics. Of course, the 
problem of stabilizing moduli and supersymmetry breaking must be addressed. However, 
it is clear that if one cannot find an MSSM-like model with this caveat, then this class 



2 



of theories can be ignoredo Finally, the cosmological constant problem would still need 
to be addressed. But perhaps the only role of the 10^'"' vacua is to resolve this problem. 

In a previous paper, "Mini-Landscape" [ML] [40] we advocated this landscape philos- 
ophy. The present paper extends the previous search and also addresses some important 
phenomenological issues. We base our model scan on the heterotic Eg x Eg string [41,42] 
compactified on an orbifold [43-49] . Our study is motivated by recent work on an orbifold 
GUT interpretation of heterotic string models [50-52]. We focus on the Zg-II = ^3X^2 
orbifold, which is described in detail in [50,52,53]. The search strategy is based on the 
concept of "local GUTs" [53-57] which inherits certain features of standard grand uni- 
fication [58-61]. Local GUTs are specific to certain points in the compact space, while 
the 4D gauge symmetry is that of the SM. If matter fields are localized at such points, 
they form a complete GUT representation. This applies, in particular, to a 16-plet of 
a local SO(IO), which comprises one generation of the SM matter plus a right-handed 
neutrino [60,61], 

16 = (3, 2)1/6 + (3, l)-2/3 + (3, 1)1/3 + (1, 2)_i/2 + (1, 1)1 + (1, 1)0 , (1.1) 

where representations with respect to SU(3)c x SU(2)l are shown in parentheses and the 
subscript denotes hypercharge (with electric charge given by Q = T3L + Y) . On the other 
hand, bulk fields are partially projected out and form incomplete GUT multiplets. This 
offers an intuitive explanation for the observed multiplet structure of the SM [53-56]. 
This framework is consistent with MSSM gauge coupling unification as long as the SM 
gauge group is embedded in a simple local GUT Giocai 5 SU(5), which leads to the 
standard hypercharge normalization I§ 

We find that the above search strategy, as opposed to a random scan, is successful 
and a considerable fraction of the models with SO(IO) and Eg local GUT structures 
pass our criteria. Out of about 3 x 10'^ inequivalent models which involve 2 Wilson lines, 
0(200) are phenomenologically attractive and can serve as an ultraviolet completion of 
the MSSM. In the present paper we extend our previous analysis in several ways. 

• In ML [40], at the last step in our analysis of a theory, we evaluated the effective 
mass operators, for the vector-like exotics, up to order 8 in fields. If all the exotics 
obtained mass, the model was retained. When calculating the rank of the mass 
matrices, we assumed that requiring the singlet configuration to respect supersym- 
metry would not change the result. 

^Note, a handful of heterotic string models with MSSM-like structure have been discussed in the 
literature [22-39]. 

^Note even if one relaxes this constraint as a prior it was shown that 90% of the MSSM-like models 
satisfying sin^ 6w =3/8 at the string scale necessarily satisfy this constraint [62]. Also, the discrepancy 
between the string scale, 0(10^^ GeV), and the 4D GUT scale, O{10^^ GeV), can in principle be resolved 
by threshold corrections due to states near the string scale. 
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In this paper we explicitly demonstrate that the decoupling of the exotics is con- 
sistent with supersymmetry. We first find the D = fiat directions. If the exotics 
decouple along these directions, then in particular models we check for F = 0. 
Then complexified gauge transformations allow us to satisfy F = and D = 
simultaneously. 

• In ML [40], we presented a model allowing for R parity. However, we did not 
perform a systematic search for R parity invariant vacua. Dangerous R parity 
violating dimension four operators can be forbidden by family reflection symmetry 
(FRS) [or matter parity] [63], i.e. a discrete Z2 subgroup of U(1)b_l (baryon minus 
lepton number). For other approaches see [21,64]. 

In this paper we evaluate B — L, searching for a "suitable" definition which has 
the accepted value on all standard model particles and gives most standard model 
singlets a value satisfying 3(i? — L) = mod 2. This condition preserves a 
subgroup of B — L under which chiral matter superfields are odd and Higgs su- 
perfields are even. Singlets with 3{B — L) = mod 2 can obtain vacuum expec- 
tation values (VEVs) for decoupling exotics, as well as giving effective quark and 
lepton Yukawa couplings. Some of these singlets give Major ana masses to right- 
handed neutrinos [65] (for earlier work see [19]), preserving R-parity. Note, if sin- 
glets with 3{B — L) = 1 mod 2 obtain VEVs, R-parity is broken and dimension 4 
baryon/lepton number violating operators are typically generated. In an Appendix, 
we also consider a possible Z^r generalization of FRS. 

The paper is organized as follows. In Section [2] we review the search strategy defined 
in ML [40]. In Section [3] we present our results solely on the issue of decoupling of 
vector-like exotics along D-flat directions. At this point we can compare our results to 
other MSSM searches in different regions of the string landscape. We show that we are 
extremely successful in finding models which have the characteristics of the MSSM. In 
the following sections we consider many of the phenomenological issues one must face on 
the road to the MSSM. In particular, in Section [4. II we discuss the problem of obtaining 
one pair of light Higgs doublets, a heavy top and then the additional constraint for a 
conserved K-pantj / family reflection symmetry. In SectionOwe discuss two models which 
satisfy the aforementioned constraints in detail. In particular we consider the effective 
Yukawa couplings for quarks and leptons in the limit that exotics decouple. We also 
study the See-Saw mechanism in these examples. Finally, in Section [6] we summarize our 
results and discuss some remaining issues. 
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2 "Mini-Landscape" search strategy [40]: Local GUTs 



Our model search is carried out in the Zg-II orbifold compactification of the heterotic 
Eg X Eg string with the twist vector 

V = ^(1,2,-3) (2.1) 
b 

acting on an G2 x SU(3) x S0(4) torus (see Fig. [U for details see [52,53]). 




Figure 1: G2 x SU(3) x S0(4) torus lattice of a Zg-II orbifold. Possible 
Wilson lines are denoted by W3, W2 and W^- The fixed points in the 
figure are those of the Ti twisted sector. 

It is well known that with a suitable choice of Wilson lines it is not difficult to obtain 
the SM gauge group up to U(l) factors. The real challenge is to get the correct matter 
spectrum and the GUT hypercharge normalization. To this end, we base our strategy on 
the concept of local GUTs. An orbifold model is defined by the orbifold twist, the torus 
lattice and the gauge embedding of the orbifold action, i.e. the gauge shift V and the 
Wilson lines Wn- We consider only the gauge shifts V which allow for a local SO(IO) or 
Eg structure, i.e. including 16— or 27— plets in the Ti twisted sector. For the twist (12. ip . 
the SO(IO) shifts are given by [66] 

yso{io),i ^ 1^ 1^ 0, 0, 0, 0, 0) (|, 0, 0, 0, 0, 0, 0, 0) , 

ySO{lO),2 ^ 1^ 1^ 0) 1^ ^ ^2.2) 

while the Eg shifts read 

y^''^ = (i h i 0' 0' 0' 0' 0) (0, 0, o, o, o, o, o, o) , 

= (|, i, i, 0, 0, 0, 0, 0) (i, i, 0, 0, 0, 0, o, o) . (2.3) 

These gauge shifts are such that the left-moving momenta P (we use the standard 
notation, for details see e.g. [51-53]) satisfying 

P • U = mod 1 , P2 ^ 2 , P G AesxEs (2.4) 

are roots of SO(IO) or Eg (up to extra group factors). In fact, this defines the "local" 
gauge symmetry in the Ti sector, for states residing at the origin in the G2 and SU(3) tori 
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and at the two fixed points in the S0(4) torus which are unaffected by the W2 Wilson 
hne along the eg direction (see Fig.[T])o The massless states of the first twisted sector are 
guaranteed to contain 16-plets of SO(IO) at the fixed points with SO(IO) symmetry or 
27-plets of Eg at the fixed points with Eg symmetry. This is established by considering 
the mass operator for left-movers, 

1 1 ^ 

_ (P + y)2 _ 1 + _ ^ 1^^1(1 _ |^.|) = , (2.5) 

i=l 

with the shift vectors V (Eqns. (12. 2p and (12. 3p ) and the twist vector v (Eq. ()2.ip ). For 
example, in the case of the two SO(IO) shifts, the massless SO(IO) spinor is given by 

P = i(-l,-l,-l,±l,±l,±l,±l,±l)(0«) (2.6) 

with an even number of minus signs. 

Since these massless states are automatically invariant under the orbifold action, 
they all survive in 4D and appear as complete GUT multiplets. In the case of SO(IO), 
that gives two complete SM generations, while in the case of Eg we have two 27s with 
27 = 16 + 10 + 1 under SO(IO). It is thus necessary to decouple all (or part) of the 
10s from the low energy theory. The third generation has to come from other twisted or 
untwisted sectors. The localized 16- and 27-plets are true GUT multiplets, whereas the 
third or "bulk" generation only has the SM quantum numbers of an additional 16-plet. 

The Wilson lines are chosen such that the standard model gauge group is embedded 
into the local GUT as 

GsM C SU(5) C SO(IO) or Eg . (2.7) 

Moreover, hypercharge is that of standard GUTs and thus consistent with gauge coupling 
unification. The spectrum has certain features of traditional 4D GUTs, e.g. matter fields 
form complete GUT representations, yet there are important differences. In particular, 
interactions generally break GUT relations since different local GUTs are supported at 
different fixed points. Also, gauge coupling unification is due to the fact that the lOD 
(not 4D) theory is described by a single coupling. 

Let us now recall the search strategy and results from the "Mini-Landscape" search 
[40]. Consider, for example, models with the SO(IO) local structure. For each of the 
SO(IO) shifts of Eq. (USD, we follow the steps: 

® Generate Wilson lines Ws and W2- 

® Identify "inequivalent" models. 

® Select models with GgM C SU(5) C SO(IO). 

^We assume that there are only two Wilson lines, W3 in the SU(3) torus and W2 in the SO(4) torus. 
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® Select models with three net (3,2). 

® Select models with non-anomalous U(l)y C SU(5). 

® Select models with net 3 SM families + Higgses + vector-like. 

The results are presented in table [TJ The models with the chiral MSSM matter content 
are listed in [67] o 



criterion 


ySO(10),l 


ySO{10),2 




yE6,2 


® inequivalent models with 2 Wilson lines 


22, 000 


7,800 


680 


1,700 


® SM gauge group C SU(5) C SO(IO) (or Ee) 


3563 


1163 


27 


63 


® 3 net (3,2) 


1170 


492 


3 


32 


® non-anomalous U(l)y C SU(5) 


528 


234 


3 


22 


® spectrum = 3 generations + vector-like 


128 


90 


3 


2 



Table 1: Statistics of Kg-II orbifolds based on the shifts 
ySO(io)a^ySO(io).2^^E6,i^yE6.2 ^j^-j^ Wilson hnes. 



To show that the decoupling of exotics is consistent with string selection rules is 
a technically involved and time consuming issue. We must select models in which the 
mass matrices for the exotics have a maximal rank such that no exotic states appear at 
low energies. We consider superpotential couplings up to order 6 in SM singlets. In our 
previous analysis, ML, we allowed any SM singlet to obtain a non- vanishing VEV. In the 
following section we refine our search and demand that all singlet VEVs be along D-flat 
directions. This requires solving the non-trivial D-flatness conditions. In this analysis we 
focus on the two SO(IO) shifts. Note, there are 218 models in this sector after step ®. 
In the following section we consider the decoupling of exotics. We do this in two steps. 
In the first step we construct the effective mass operators for the exotics and check to 
see if the exotics decouple allowing arbitrary singlet VEVs. In the second step we only 
consider singlet VEVs along D-flat directions. 

^In a recent paper [62], the two constraints, (D SM gauge group C SU(5) C SO(IO) (or Ee) and © 
non-anomalous U(l)y C SU(5) were removed. This search has lead to about 10 times more models. 
However the additional constraint that sin'^ 9w = 3/8 reduced this number by 90% so that there were 
only a handful of additional models. It suggests that in order to find the MSSM, one may need to require 
local GUTs. 
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3 Decoupling exotics 



We evaluate all effective mass operators for the exotics to order 6 in SM singlet 

fields 3j, 

W D Xj (si • • -Jjv) ■ (3.1) 

In general, s transform non-trivially under the extra U(l)s and hidden sector gauge 
groups. To construct the mass operators ()3.ip . we find all monomials of the above form 
consistent with string selection rules. These rules have been discussed previously in the 
literature. They include space group and R-charge selection rules, in addition to the 
standard field theoretic requirement of gauge invariance. Complete details of these string 
selection rules are given in Appendix A. We should emphasize here that in previous 
analyses a 7 selection rule has also been enforced [52,68,69]. We disagree with this 
additional 7 rule and in Appendix A we give a general argument why this rule is not a 
selection rule§| 

We consider the 218 models remaining after step ® from the two SO(IO) shifts (128 
from ySO{iO),i 90 fj-om ySO{lO),2^ 

see Table [TJ If in a particular model all exotics 
decouple to order 6 in the product of s fields, assuming arbitrary s VEVs |§, we retain 
the model. The number of models satisfying decoupling at this step is 191 (106 from 

ySO{10),l j^^j 85 f^Qj^ ySO{10),2)_ 

now determine D-flat directions for all s fields. 
Our procedure for determining D-flat directions is described in Appendix B. We then 
retain the subset of the 191 models for which the exotics decouple along D-flat directions 
to order 6 in the s fields. We find 190 models remaining. Clearly, D-flatness does not 
impose an important constraint. Thus we are successful in 190/3 • 10"^ or 0.6% of the 
cases. 

The results of our search may now be compared to many other searches in the 
literature. We have 218 models with the SM gauge group, 3 families and only vector- 
like exotics from our two SO(IO) shifts. Out of these we find 190 for which all exotics 
decouple along D-flat directions. In certain types of intersecting D-brane models, it was 
found that the probability of obtaining the SM gauge group and three generations of 
quarks and leptons, while allowing for chiral exotics, is less than 10~^ [8,9]. The criterion 
which comes closest to the requirements imposed in [8, 9] is ®. We find that within our 
sample the corresponding probability is 6 %. In [6,7], orientifolds of Gepner models were 
scanned for chiral MSSM matter spectra, and it was found that the fraction of such 
models is 4 x 10"^'^. These constructions contain the MSSM matter spectrum plus, in 
general, vector- like exotics. This is most similar to step ® in our analysis where we find 
218 models out of a total of 3 x lO'^ or 0.7 %. In comparison, approximately 0.6 % of our 

^In fact, we have shown that all exotics decouple in model Al in Ref. [52] if one eliminates the 7 rule. 
®Note that giving VEVs to the s fields can often be interpreted as blowing up the orbifold singularities 
(for recent developments in this direction see [70-73]). 
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models have the MSSM spectrum at low energies with all vector-like exotics decoupling 
(with exotic mass terms evaluated to order s^) along D-fiat directions. Note also that, 
in all of our models, hypercharge is normalized as in standard GUTs and thus consistent 
with gauge coupling unification. 

4 Road to the MSSM 

In this section we consider other phenomenological hurdles which must be overcome in 
order to reach the MSSM. These hurdles include finding supers3mimetric minima with 
proton stability, an exactly conserved R-parity, a fi term for the light Higgs doublets of 
order the weak scale, a top quark Yukawa coupling of order 1, gauge coupling unification, 
and more. 

4.1 Constraints 
R-parity conservation 

One of the most formidable obstacles in string constructions is obtaining a conserved 
R-parity. In this paper we propose one possible route, i.e. obtaining a "family reflec- 
tion symmetry" or "matter parity". In this regard, wc evaluate B — L, searching for 
a "suitable" definition which has the accepted value on all standard model particles, 
is vector-like on all exotics and produces a number of SM singlets with even and zero 
3{B - L) charge. 

Giving such singlets VEVs preserves a subgroup of 5 — L, denoted family reflec- 
tion symmetry or matter parity, under which chiral matter superfields are odd and Higgs 
superfields are even. We find that the exotics can be decoupled and the right-handed 
neutrinos can be given Majorana masses consistent with this symmetry. In Appendix C, 
we show that it is possible to allow any s field to obtain a VEV as long as it has B — L 
eigenvalue / = 0, ±2/ (2Z -|- 1). This will leave invariant 1,2^. 

To apply the above strategy, we must first give a "suitable" definition of -B — L. A 
possible algorithm to identify the corresponding generators is discussed in Appendix D. 
Upon defining B — L, we must verify D-flatncss for the subset of SM singlets with B — L 
charges / = 0, ib2/(2Z + 1) and check that all exotics decouple. This is a tedious task, 
requiring much computer time. In order to minimize the amount of time, we focus our 
attention on a subset of the 190 models which have renormalizable top quark Yukawa 
couplings. 
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Light Higgs doublets 



The Higgs doublets of the MSSM are vector-hke and genericahy in our analysis all the 
Higgs doublets decouple. Retaining one pair of light Higgs doublets in the MSSM is the 
/i problem, and we must now face this issue. We look for vacuum configurations in which 
the fi term vanishes to a certain order in the s fields. At the same time we require that 
all the exotics decouple. Of course, it would be nice to have a symmetry argument for a 
small /Lt. 

Order one top quark Yukawa coupling 

The top quark Yukawa coupling is necessarily of order 1. Hence it is natural (although 
perhaps not absolutely necessary) to require that for the top quark we have a renormal- 
izable 0(1) Yukawa coupling (3, 2)i/g (3, l)_2/3 (1, ^)i/2: i-^- one of the following types 

UUU, UTT , TTT , (4.1) 

where U and T denote generic untwisted and twisted fields, respectively. The UUU cou- 
pling is given by the gauge coupling, UTT is a local coupling and thus is unsuppressed, 
while the TTT coupling is significant only when the twisted fields are localized at the 
same fixed point. We discard models in which the above couplings are absent or sup- 
pressed. In ML we required that the top quark have a Yukawa coupling at tree level, i.e. 
cubic order in the fields, in addition to decoupling of all exotics, albeit assuming arbitrary 
VEVs for the SM singlets. Of the 190 models which decouple along D-flat directions we 
have 105 (85) coming from the first (second) SO(IO) shift. Out of these there are 55 
(32) with "heavy top" and 50 (53) with "no heavy top." We thus find 87 models which 
decouple along L'-flat directions and have a "heavy top." Note, this is just one less than 
discussed in ML at step ®. 

R-parity invariant models with cubic top Yukawa coupling 

We find a "suitable" definition oi B — L for 34 of the 55 (5 of the 32) models of the 
first (second) SO(IO) shifts. Note however that for each case there are several possible 
inequivalent choices. This is because of two ambiguities which need to be resolved. 

1. In many cases there are vector- like exotics with SM gauge charges identical to those 
of quarks, leptons and Higgs doublets. Thus there arc difi"crcnt ways to choose which 
of these states have standard B — L charges. Each choice can lead to a different 
definition oi B — L. 

2. For each choice of SM particles above, there may be more than one B — L definition. 
In some cases there are continuous families of solutions. 
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Including all of these possibilities we find 3447 (144) suitable B — L generators from 
the first (second) SO(IO) shifts, which also lead to the presence of SM singlets with 
charges B - L = 0, ±2, ±2/3, ±2/5, . . . ,±6/7. We find, however, that these lead to 85 
(8) inequivalent models. Requiring the absence of extra unbroken U(l)s reduces this set 
to 42 (0) acceptable models. Finally, demanding that all exotics decouple along D-&at 
directions leads to 15 (0) acceptable solutions with an exact low energy R-parity. This 
result is specific to our (S — L)-based strategy and we expect, in general, more acceptable 
models to exist. 

4.2 Approaching the MSSM 

Further issues to be addressed are as follows. 

1. We must check that quarks and leptons obtain non-trivial masses. For neutrinos, 
this includes an analysis of Majorana masses and the See-Saw mechanism. 

2. We must also consider dimension 5 baryon and lepton number violating opera- 
tors [63,74,75]. These operators are not forbidden by R-parity and are typically 
generated. Their coefficients must necessarily be suppressed in order to be consis- 
tent with proton decay experiments [76,77]. 

3. Precision gauge coupling unification should be addressed [20,52,78]. This includes 
a calculation of the string threshold corrections [80-82] . 

4. Finally, F = has to be verified. This constraint guarantees that our vacua are 
indeed supersymmetric. In general, F = solutions exist. Some of them can be 
found numerically by truncating the superpotential and solving polynomial equa- 
tions. Once they are found, F = and D = can be satisfied simultaneously using 
complexified gauge transformations [83] (for a detailed discussion see [53]). 

All of these checks are clearly time consuming and we have not performed an inclusive 
analysis. We have however found many vacua with R-parity. In the next section we discuss 
our results for two particular examples. In these examples we have demanded that: 

• all exotics are massive, 

• there is one pair of massless Higgses, 

• the mass matrix for the right-handed neutrinos has full rank, 

• no extra U(l) factors remain, 

• hidden sector gaugino condensation is possible, 

• R-parity is unbroken. 
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5 Two "Benchmark" models with R-parity 



We now discuss two particular "benchmark" models. In Model 1 we also consider two dif- 
ferent vacuum configurations and show how the phenomenology depends on the different 
choices of vacua. 

5.1 Model 1 

The model is defined by the shifts and Wilson lines given in Appendix E. The gauge 
group after compactification is 

Gorbifoid = SU(3) X SU(2) X [SU(4) x SU(2)'] x \J{lf . (5.1) 

The resulting massless spectrum includes three SM generations plus vector-like ex- 
otics with respect to the SM gauge group. 

The model allows us to define a "suitable" B — L generator, 
tB-L = (^0,0,0,0,0, (0,0,0,0,0,2,0,0) . (5.2) 
with two essential properties (cf . Table [2|) : 

• the spectrum includes 3 generations of quarks and leptons plus vector-like exotics 
with respect to Gsm x U(1)b-l , and 

• there are SM singlets with B — L charge ±2. 

In the following discussion we consider two different vacuum configurations to illus- 
trate the dependence on the particular vacuum class, i.e. the set of SM singlets with 
non-zero supersymmetric VEVs. 

5.1.1 Model 1, vacuum configuration A 

Consider a vacuum configuration where the fields 

{Si} = {Xl, X2, X3, X4, /ll, /l2, /is, hi, /l5, /l6, /19, /llO, ^S' 4^ ^IH '^IS' '^IS' 

„0 „0 „0 „0 „0 „0 „0 „0 „0 „0 „0 „0 T. /-c: o\ 

*16' *17' '^18' '^20' *21' *22; *23' '^25' *26' '^27' *30' *31J K'^-'^) 

develop a VEV while the expectation values of all other fields vanish. In this vacuum 
configuration we set 14 of the original 46 SM and hidden SU(4) singlets to zero. The 
emerging effective theory has the following properties: 

1. the unbroken gauge symmetries are 

Gsm X Ghid , (5.4) 
where Ghid = SU(4). 
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Table 2: Spectrum. The quantum numbers under SU(3) x SU(2) x [SU(4) x 
SU(2)'] are shown in boldface; hypercharge and B — L charge appear as 
subscripts. Note that the states sf, fi, fi and mi have different B — L 
charges for different i, which we do not explicitly list. 



2. since B — L is broken by two units, there is an effective matter parity Z2 

3. there is only one pair of Higgs candidates, (pi and 0i; the /i-term 



M 



(5.5) 



<f>i=(pi=0 



vanishes up to order 3*, at which we work. That is, there is one pair of massless 
Higgs doublets. 

4. we check that the solution satisfies F = for all fields (cf. F-flatness below) and in 
addition, switching on {3j}-fields allows us to cancel the FI term without inducing 
D-terms (cf. Appendix iBl). 



5. all vector-like exotics decouple (cf. Appendix IE. 2p . 

6. neutrino masses are suppressed via the see-saw mechanism. 

That is, we have obtained a super symmetric vacuum with the precise matter content 
of the MSSM with an exact R parity. This has to be contrasted to [53] where R parity 
was approximate, and to [33] where R parity exists only at the classical level, and where 
yj= and Y(i vanish at the same level. Our model also does not suffer from the problem 
encountered in [37], where it was found that you can either decouple all exotics or have 
i?-parity but never both. 
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Charged fermion Yukawa matrices 

The charged fermion Yukawa matrices ar^ 



Y„, 



f 

^ ^ ^ \ , Yh 

V 56 S6 1 



/ 















I 








Y. 








~6 


35 








36 


~6 






(5.6) 



Here, for n a non-negative integer, represents the the smallest value of n for which this 
term appears in the matrix; thus giving the dominant contribution. Each term is in fact 
a sum of monomials containing several different SM and hidden SU(4) singlets. The up- 
type quark Yukawa matrix is given directly in terms of the coupling of the up-type Higgs 
to the three q and u fields. The down-type quark and charged lepton Yukawa matrices 
are obtained by integrating out a pair of vector-like d- and d-quarks and i- and ^-fields, 
respectively. We find that the up and charged lepton Yukawa matrices have rank 3, while 
the down quark Yukawa matrix has only rank 2, at this order in s singlets. However, we 
have checked that at order 8 in 3 fields has rank 3. 

Neutrino matrices 

In our vacua SU(2)' is broken such that the SU(2)' doublets iji and fji correspond to SM 
singlets with qb-l = ilj 







m = { _ 1 = ( _ I (5.7a) 

and 

Vi = 

The dimensions of the "right-handed" neutrino mass matrices are 

Mnn = 23 X 23 , 

Mnn = 23 X 26 , 

Mnn = 26 X 26 , (5.8) 



with the complete v — ly mass matrix given by 

Mnn M 

M^- M 



^Quark and lepton doublets multiply the Yukawa matrices on the left. 
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We have checked that it indeed has full rank. For more details, see Webpage [67]. 
The Dirac neutrino Yukawa couplings have dimensions 



n 



= 4 X 23 



Yn = 4 X 26 . 



(5.10) 



The effective light neutrino mass operator emerges from 



(5.11) 



where Yi, = {Yn, Yn), by integrating out the pair of heavy leptons £i and i' where £' is a 
linear combination of the ii . We have checked that the light neutrinos all obtain a small 
mass. The large dimension of the matrices effectively reduces the See-Saw scale [65]. 
Also, neutrino phenomenology works differently in the presence of many 'right-handed' 
neutrinos [84,85]. 

Dimension 5 baryon and lepton number violating operators 

We further analyzed the question of dimension 5 proton decay operators. We find that 
both qqql and uude appear at order S^. They are also generated by integrating out the 
heavy exotics. For example, the following couplings exist 

qi h 64 , qi h 65 , q2 h ^4 , 92 h h , qi qi , qi qi 65 , q2 q2 , q2 S5 ■ (5.12) 

Hence integrating out the states 6i, 6i produces dangerous dimension 5 operators. These 
must be sufficiently suppressed to be consistent with present bounds on proton decay 
[76,77]. We have verified that, for some particular sVEVs, it is possible to suppress the 
qqqi operators induced by the trilinear couplings (j5.12p . However, higher order couplings 
also introduce baryon and lepton number violating operators. We have not been able to 
identify a suppression mechanism for such operators yet. 

/x-term and Minkowski space 

Since our singlet configuration satisfies F = D = 0, the vacuum energy is zero in the 
global SUSY limit. In supergravity, one should include non-perturbative moduli poten- 
tials which would be responsible for spontaneous SUSY breaking. In fact, in the context 
of gaugino condensation [86-89], the SU(4) subgroup of the second Eg gives rise to TeV 
soft masses, which is a common feature of our "fertile" patch of the landscape [95]. A 
specific realization of SUSY breaking via gaugino condensation in heterotic string com- 
pactifications is given by Kahler stabilization [90-92]. Requiring a Minkowski vacuum 
puts a constraint on the total superpotential which includes contributions from charged 
matter and moduli. 
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An interesting feature of our singlet configuration is that the charged matter super- 
potential vanishes at order 6 in singlet fields, 



{W{s}) = 0. (5.13) 

In fact the superpotential is a polynomial in s fields and in this vacuum configuration 
each monomial term in W vanishes independently. Therefore, the total superpotential is 
given solely by its non-perturbative part. This is expected to be very small and thus a 
small gravitino mass and a small cosmological constant can in principle be achieved. 

In this model we also find an intriguing correlation between the /x-term and W^s). 
Since the Higgs doublets are untwisted and the combination (j)i has vacuum quantum 
numbers, requiring that each monomial appearing in ^ vanishes also implies (W(s)) = o|f| 
This means that the /U-term is of the order of the expectation value of W , i.e. the gravitino 
mass. 

F-flatness 

We analyze the F-terms in the configuration defined by (j5.3p . The only non- vanishing 
F-terms are 

= % ^^^"^^ ^ i"3> S%, S?4, s%} . (5.14) 
They read 

F3 = c,27 4^4+a2ss%(s0,)^s0+a29s0M^sl, (5.15a) 
Fr = 041 hihios^gs1+a44, h2hgs'^gs1+a4e, s^^s'^c^s'^Qsl+air s^s^^s'^gs'^gsl+aso *i8*27''30*4+"55 «i5(^*3o)^^4 
+"56 46i4of4+'^57 S^oS%sO^S°+a5S ''^i4o4l'^2 ■ (5.15b) 
Fio = ai hihios'lgS%+a4 h2hQS%s%+ai2 s%s1^s'^^s%+ais s1s1^s%s^gS%+aie Si3«i8^274o 
+021 S%S°^{s%)^+a22 Sl3'5?6(4o)^+"23 Sl3«2o4o4l+"24 ^IS^Sl^O'^Sl+^Sl ^26^9 
+as2 sU4e)'4+'^S3 s%{s°,,fsO , (5.15c) 
Fi4 = agg /ll?tl0S30«9+a92 h2h<)s'^gs'^+ag4 sl^s'i^^s'^gS^+ags S?S?gS26^30''9+"98 «18'527'530'5g 

+«103 *?5(4o)^'^9+«104 S?6(4o)'^9+"105 '^2o4o4l''9+°106 ^^M^O^l^ . (5.15d) 
Fig = 035 /ll?t2S3gSgS^+a36 s5-,s5gs!^gS^S^+a37 S2qS22*30*5'^6~'"'-'38 ■^21'^22'^30*5*6"'"'-'39 *20^23'^30'^5*6 

+^40 s'^^sO^s^sOsO. (5.15e) 

Here a, denote superpotential coefficients. The F-term equations, Fi = 0, have trivial 
and non-trivial solutions. For instance, (|5.15ap has the trivial solutions 

4 = or si = 0, 



This applies to the untwisted Higgs pairs in many models of our MiniLandscape, for instance also to 
the model presented in [53,55]. 
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as well as the non-trivial solution relating various VEVs, 

"29 ^26 

The strategy is now to take the non-trivial solution and insert it into the other equations. 
By doing so, one can trade the i^-term constraints for relations between the VEVs. We 
find that this strategy is successful and we can satisfy all equations with non-trivial s 
VEVs. 



5.1.2 Model 1, vacuum configuration B 

Now consider the 'vacuum' configuration where the fields 

{si} = {xi,X2, X3, X4, hi,h2, h3, hi, h^, he, hj, hs, kg, hw, s?, s^, s^, si, 

„o „o „o „o „o „o „o „o „o „o „o „o „o „o „o 

*5' *6' *8' '^9' *10' *12; *13' *14' '^15' '^16' *17' '^18' *20' 

„0 „0 „0 „0 „0 „0 „0 „0 „0 „0 „0 „0 T. f^-\7\ 
*2l5 *22; '^23' '^24' '^25' '^26' *27' *28' *29' '^30' *31' *32J ) 

develop a VEV while the expectation values of all other fields vanish. Hence, in this 
vacuum configuration only one SM singlet VEV is set to zero, i.e. (s^g) = 0. The 
emerging effective theory has most properties identical to those in Model lA: 

1. the unbroken gauge symmetries are 

GsM X Ghid , (5.18) 
where Ghid = SU(4). 

2. since B — L is broken by two units, there is an effective matter parity Zg^. 

3. there is only one pair of Higgs candidates, (pi and 0i; the /i-term 



(5.19) 

4>i=<t>i=o 



vanishes up to order s*, at which we work. That is, there is one pair of massless 
Higgs doublets. 

4. we check that the solution satisfies F = Q for all fields and in addition, switching 
on {si}-fields allows us to cancel the FI term without inducing D-terms. 

5. all vector-like exotics decouple (cf. Appendix IE. 3^ . 

6. neutrino masses are suppressed via the see-saw mechanism. 

However, the detailed form of the Yukawa and exotic mass matrices has changed. 
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Charged fermion Yukawa matrices 



The charged fermion Yukawa matrices are 



Y, 



J5 J5 J5 



\ 


0/ 



Y. 




(5.20) 



The up-type quark Yukawa matrix is given directly in terms of the coupling of the up- 
type Higgs to the three q and u fields. The down-type quark and charged lepton Yukawa 
matrices are obtained by integrating out a pair of vector-like d- and d-quarks and i- 
and ^-fields, respectively. We find that (just as in Model lA) the up and charged lepton 
Yukawa matrices have rank 3, while the down quark Yukawa matrix has only rank 2, at 
this order in s singlets. In fact, to this order in SM singlet fields, the superpotential does 
not couple two right-handed down quarks, ^3^4, to the quark doublets. This is because 
c^3,4 are in the T4 twisted sector. However, we have verified that some of the zeros in 
Yd get filled in at higher orders and at order 8 has rank 3. Note that in this vacuum 
configuration the Yukawa matrices retain a form consistent with the underlying D4 family 
symmetryl§ 



Neutrino matrices 



We have checked that all right-handed neutrinos obtain mass in this vacuum configura- 
tion. Thus the See-Saw mechanism works exactly as in Model lA, although the detailed 
forms of the matrices differ. For more details, see Webpage [67]. 



Dimension 5 baryon and lepton number violating operators 

We further analyzed the question of dimension 5 proton decay operators. We find that 
both qqqi and uude appear at order s^. They are also generated by integrating out the 
heavy exotics. For example, the following couplings exist 

qi ii ^4 , qi li 4 , q2 ^2 (^4 , 92 ^2 (^5 , Qi Qi ^4 , Qi Qi S5 , q2 Q2 , q2 q2 S5 ■ (5.21) 

Hence integrating out the states 6i, 6i produces dangerous dimension 5 operators. These 
must be sufficiently suppressed to be consistent with present bounds on proton decay 
[76,77]. We have verified that, for some particular s VEVs, it is possible to suppress the 
qqqi operators induced by the trilinear couplings (|5.2ip . However, higher order couplings 

^The D4 family symmetry is a consequence of the space group selection rules and the geometry of the 
SO(4) torus (Fig. [l| [52,93]. States sitting at the two vertical fixed points on the S0(4) torus transform 
as doublets under D4- 
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also introduce baryon and lepton number violating operators. We have not been able to 
identify a suppression mechanism for such operators yet. 



^-term and Minkowski space 

As in Model lA, requiring a Minkowski vacuum puts a constraint on the total super- 
potential which includes contributions from charged matter and moduli. An interesting 
feature of the present vacuum configuration is that the SM matter singlet superpotential 
to order is of the form, 

W(S) = J2Pi{D)Pi(S), (5.22) 

i 

where Pi are polynomials in SM singlet fields (the index i labels a particular polynomial) 
which are either D4 doublets, which we now re-label as -D, or SM and D4 singlets, s. In 
particular, the D4 doublets which enter W(s) are 

D, = {slsl) D2 = {sl,s%) (5.23) 

55 = (s?,s?3) De = isl,s^,^) . 

The polynomial in D4 doublets is, to this order, quadratic in doublets and is given by 
the trivial D4 singlet scalar product, for example, 

Di-D2 = (si si + si s%) . (5.24) 

We then find (up to calculable dimensionful coefficients in units of the string scale) 



+ (^Di ■ De + D2 ■ 4 



30 *3o(*15 + ^le) + ^17(^25 + *28) 

+ •''■l8('^24 + '^27) + '53l('520 + ^21) + ^32(^22 + '^23) 

+ {S% + S$S$8 + S^s27)(46 + S29) + ^1(^8 + hio) + h2{h7 + hg) 

+ • D4) S%slo (s°75?8 + hih2 + (s°o + 4l)(^22 + 4)) • (5.25) 

Thus, to order 6 in SM and hidden SU(4) singlets, the polynomials Pi{D) are completely 
determined by the D4 symmetry, while the polynomials Pi{s) arc non-trivial for all i. 
One particular F = D = solution is given by the roots of {Pi{D)) = {Pi(s)) = for all 
polynomials i. Hence, once again, 

{W{s)) = 0. (5.26) 
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Therefore, the total superpotential is given solely by its non-perturbative part. This is 
expected to be very small and thus a small gravitino mass and a small cosmological 
constant can in principle be achieved. 

In addition, just as in Model lA, the //-term contains all terms present in W(s) and 
to order 6 in SM singlets we have /i = 0, when W = 0. This means that the //-term is of 
the order of the expectation value of W, i.e. the gravitino mass. 



5.2 Model 2 

The model is defined by the shift and Wilson lines given in Appendix F. It was already 
included as an example in ML [40] and [65]. 

The gauge group after compactification is 
G = [SU(3) X SU(2)] X [S0(8) x SU(2)] x U(l)^ . (5.27) 

As before, we are able to define a suitable B — L generator, 

/ 2 2 2\ /I 1 1 1 \ , , 

ts-L = (l,l,0,0,0,--,--,--j (2'2'°'2'2'°'°'°J ^'-^^^ 

with two important properties (cf. Table [3]): 

• the spectrum includes 3 generations of quarks and leptons plus vector-like exotics 
with respect to Gsm x U(1)b-l , and 

• there are SM singlets with B — L charge ±2. 

Consider a 'vacuum' configuration where the fields 
r-r.i _ r^, ^„ ^0 ^0 „o „o „o „o „o „o „o „o „o „o „o 

\°iS — lAl; A2, ■^5, ^8' *9' '^125 *15' '^16' *22; '^24' '^35' *41> *43' '^46' 

/l2, /l3, /l5, ^9, ^13, /ll4, /l20, /i21, /i22} (5.29) 

develop a VEV while the expectation values of all other fields vanish. The emerging 
effective theory has the following properties: 

1. the unbroken gauge symmetries are 

Gsm X Ghid , (5.30) 
where Ghid = SO (8). 

2. since B — L is broken by two units, there is an effective matter parity ^2^. 



20 



it 


irrep 






44- 

it 


irrep 


Inlipl 
IdUfc,! 


3 


(3 2- 


1 1' 


(1/6,1/3) 


HI 




3 


(3 1 




) (-2/3,-1/3) 


III 


O 






(1,1) 






4 


(1 9- 

I-*-, 




(0,*) 




4 




1 1 


) (1/3,-1/3) 


J- 




1 
± 






(-1/3,1/3) 


rl- 


4 


(1,2; 


1, ly 


(-1/2,-1) 



H 




1 
1 


(1,2; 


1, 1 


'(1/2,1) 


I 
H 


4 


(1,2; 


1, 1; 


(-1/2,0) 


(pi 




4 


(1,2; 


1 1 

1, 1 


) (1/2,0) 


(pi 


q 
o 


(«>,l 


1 1 

1, 1 


) (1/3,2/3) 






q 



(«>, 1; 


1 1 

1, 1 


)(-l/3,-2/3) 


X 

(>i 


20 


(1,1; 


1,1; 


(1/2,*) 


+ 




20 


(1,1; 


1,1 


)(-l/2,*) 


^i 


io 


(1, 1, 


1, ly 


(0,1) 


fii 






(1, 1, 


1, 1 


'(0,-1) 


rii 


q 


(1,1; 


1,2; 


(0,1) 






q 



(1,1; 


1,2 


(0,-1) 


Vi 


20 


(1,1; 


1,2] 


(0,0) 


1 

hi 




2 


(1,2; 


1,2 


(0,0) 


Vi 


2 


(1,1; 


1,2: 


(1/2,1) 


4 




2 


(1,1; 


1,2 


(-1/2,-1) 


^7 


2 


(1,1; 


1,1' 


(0,±2) 


Xi 




18 


(1,1; 


1,1 


(0,0) 




4 


(3,1 


1,1 


)(-l/6,*) 






4 


(3,1; 


1,1 


(1/6,*) 




2 


(1,1; 


8,1; 


(0,-1/2) 


h 




2 


(1,1; 


8,1 


'(0,1/2) 




5 


(1,1; 


8,1: 


(0,0) 


Wi 















Table 3: Spectrum. The quantum numbers under SU(3) x SU(2) x [SO (8) x 
SU(2)'] are shown in boldface; hypercharge and B—L charge appear as 
subscript. Note that the states sf, rrii and Vi have different B-L charges 
for different i, which we do not explicitly list. 



3. the Higgs mass terms are 



Jij Yj 1 



where 











s 
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J 



The up- type Higgs hu is a linear combination of 01, ^3 and ^4, 

hu ~ S^(pl + 03 + 04 , 

while the down-type Higgs is composed out of 02 and 03, 

hd ~ 02 + 03 ■ 



(5.31) 



(5.32) 



(5.33) 



The vacuum configuration is chosen such that the /x-term, being defined as the 
smallest eigenvalue of A^t , , 



dhd dhu 



(5.34) 



h„=hd=0 



vanishes up to order S*, at which we work. 
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4. we check that switching on {sj}-fields allows us to cancel the FI term without 
inducing D-terms (cf. Appendix IB]) . 

5. all exotics decouple (cf. Appendix IF.2|) . 

6. neutrino masses are suppressed via the see-saw mechanism. 

Thus, again we have obtained a supersymmetric vacuum with the precise matter 
content of the MSSM and R parity. 



Charged fermion Yukawa matrices 

The up-Higgs Yukawa couplings decompose into 

4 

^Yukawa ^ '^{yu)if Qi Uj 4>k , 

where 



k=l 
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Thus, the physical 3x3 up-Higgs Yukawa matrix is 
Y. ~ i^Y«+Yi3)+j4YW = 




(5.35) 



(5.36a) 



(5.37) 



Note that due to the Higgs mixing the top quark Yukawa coupling for this vacuum 
configuration is given by s^. Thus the corresponding s VEVs are required to be quite 
large. 

The down-Higgs Yukawa couplings decompose into 

4 



Wvukawa ^ ^(Xd)\yqid 



(5.38) 
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(5.39a) 
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(3) 
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The physical 3x3 down-Higgs Yukawa matrix emerges by integrating out a pair of 
vector-Uke d— and d— quarks, 

-S3 



/ 1 \ 

S-4 



(5.40) 



\ s / 

We note that both the up and down quarks are massless at order 6 in SM singlets. 
However, we have checked that the up quark becomes massive at order 7 and the down 
quark gets a mass at order 8. 

The charged lepton Yukawa couplings decompose into 

4 



WVukawa D l^(^e) 
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(5.41) 
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(5.42a) 
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The physical 3x3 matrix emerges by integrating out a pair of vector-like I— and 
leptons, 
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(5.43) 



Neutrino masses 

We consider vacua where SU(2)' is broken. This means that the rji and fji give rise to 
further SM singlets with qb-l = ±1, 



ni7 



^20 
"21 



and r]i 
23 



ni3 



ni7 
ni8 



(5.44) 



The dimensions of the "right-handed" neutrino mass matrices are 



Mnn = 18 X 18 , (5.45) 
Mnn = 18 X 21 , (5.46) 
Mnn = 21x21, (5.47) 

with tlie neutrino mass matrix given by 

\ Mnn Mnn J 

We have checked that it has full rank. 

The neutrino Yukawa couplings decompose into 

4 

^Yukawa D ^(i^^)!,^^ 4 n.^jfc + (F^,)g) 4 n^^fe, (5.49) 
k=l 

where yJi^\y^^^ are non- vanishing 4 x 18 matrices, yJ^''^^^ = and Yi^\YP,YP,YP 
are non-vanishing 4 x 21 matrices. The effective neutrino mass matrix obtained as 

K = Y^M^^yJ , (5.50) 

where Y,y = {Yn,Yn) and k has non-zero determinant. See Webpage [67] for details. 

Dimension 5 baryon and lepton number violating operators 

We have looked for effective dimension 5 baryon and lepton number violating operators 
in this model. We find that to order no such operators exist. However, these operators 
can be generated once the exotics 5j, 5j are integrated out. Fortunately, a clever choice 
of VEVs for the fields {si} can guarantee sufficient suppression of all induced qqqi 
operators, consistent with current bounds on proton decay [76,77]. 



^-term and Minkowski space 

Unlike in the previous model, there is no relation between the ^-term and W{s). This 
is because the Higgs doublets do not come entirely from the untwisted sector. Requiring 
spontaneous SUSY breaking in a Minkowski vacuum puts a constraint on the moduli 
VEVs. Fine-tuning is likely to be necessary to obtain a realistic gravitino mass as well 
as a small cosmological constant. 
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6 Conclusions and Discussion 



In this paper we have described the construction of heterotic MSSMs with R parity. Our 
setup is based on a particular Zg-II orbifold with an SO(IO) local GUT structure. In 
the first part of the paper we have obtained 218 models with the MSSM gauge group 
structure, 3 light families and vector-like exotics. We show that all the vector-like exotics 
can decouple along L>-flat directions for 190 of these models^ The total number of 
inequivalent models with SO(IO) shifts and 2 Wilson lines is 3 • 10^. Hence 0.6% of our 
total model set are MSSM candidates. This can be compared with D brane constructions 
where the probability of getting MSSM-like models is much less than 10~^ or Gepner 
orientifold constructions where this probability is 10~^^. 

In the second part of the paper we go further down the road towards the MSSM. 
We define a successful strategy for obtaining models with an exact R parity. We find 87 
models which have a renormalizable top Yukawa coupling. We identify 15 models with 
an exact R parity, no light exotics or U(l) gauge bosons and an order one top quark 
Yukawa coupling!^ 

We present two explicit "benchmark" examples satisfying the following criteria: 

• MSSM spectrum below the string scale - 

all exotics decouple; 

one pair of light Higgs doublets; 

top quark Yukawa coupling of order 1; 

non-trivial Yukawa matrices for charged fermions; 

See-Saw mechanism for neutrinos; 

• an exact R parity. 

The two examples have different phenomenological properties such as different struc- 
tures of the Yukawa coupling matrices and dimension 5 operators. In particular, the 
Yukawa matrices Yd and have more non-vanishing entries in Model IB than in Model 
lA. In both Models lA/B the lightest down type quark is massless at order 6 in SM 
singlets and becomes massive at order 8. The top Yukawa coupling is order one in Mod- 
els lA/B, while it is order 3^ in Model 2. This is due to Higgs doublet mixing in the 

""^"in this analysis, we have taken into account superpotential terms up to order 6 in SM singlets. At 
higher orders, we expect more models to be retained. 

^^The number 15 is a lower bound, since our search is based on a specific strategy related to B — L 
symmetry. Furthermore, more models are retained if we do not insist on having a renormalizable top 
Yukawa coupling. Also one can drop the strict constraint that exotics be vector-like with respect to B — L. 
For example, two exotics x, x with B — L charge -1 can get mass from a SM singlet VEV with charge 
+2. 
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latter. In Model 2, both the up and down quarks are massless at order 6 in SM singlets. 
However, the up quark becomes massive at order 7 and the down quark gets a mass at 
order 8. 

An interesting feature of Model 1 is that there is a correlation between the /x term 
and the expectation value of the superpotential. In fact the pair of Higgs fields are the 
only vector-like fields whose mass is correlated with the expectation value of the super- 
potential, while all exotics can consistently get mass with W = 0. This provides a novel, 
stringy solution to the MSSM /x problem. Indeed, in Models lA/B the vacuum expec- 
tation value of the superpotential and fi both vanish at order 6 in SM singlets. Thus, 
neglecting non-perturbativc effects, this model leads to a supersymmetric Minkowski vac- 
uum with 11 = 0. One expects that when non-pcrturbativc effects (hidden sector gaugino 
condensation) are taken into account, supersymmetry is broken at a hierarchically small 
scale and, because of the correlation between and {W), /x is of order the gravitino 
mass. In Model 2, on the other hand, the superpotential does not vanish in this limit 
and inclusion of non-perturbative contributions to the superpotential is necessary. 

Dimension 5 baryon and lepton number violating operators come from two sources. 
They are generated in the superpotential to some order in SM singlets. They may also 
be generated when integrating out heavy exotics. In Models lA/B the direct dimension 
5 operators appear at order s^, while in Model 2 they do not appear at this order. In 
addition, in Models lA/B and Model 2 dimension 5 operators appear when integrating 
out heavy exotics. In Model 2 these can be sufficiently suppressed with some fine-tuning. 

There are some phenomenological issues that we have not addressed in this paper. 
In particular, we have not studied precision gauge coupling unification. Although hyper- 
charge is normalized as in 4D GUTs thus allowing gauge coupling unification in the first 
approximation, there arc various corrections that can be important. First, a detailed 
analysis would require the calculation of string threshold corrections in the presence 
of discrete Wilson lines. However in specific cases these corrections are known to be 
small [94]. Second, there are corrections from the vector-like exotic states. It is possi- 
ble that precision gauge coupling unification may require anisotropic compactifications, 
leading to an effective orbifold GUT [20,52,78,79]. 

Another issue concerns proton stability. The examples we studied are challenged 
by the presence of dimension 5 proton decay operators. Their suppression may require 
additional (discrete) symmetries. There are also dimension 6 operators, generated by 
GUT gauge boson exchange, which we have not discussed. 

Finally, there are the usual questions of moduli stabilization and supersymmetry 
breakdown in a Minkowski vacuum. Some of them we discussed previously in [95]. We 
have not addressed all of these issues here. On the other hand, it is clear that if given the 
freedom of arbitrarily tuning moduli VEVs we are not able to find the MSSM, the whole 
approach would be futile. However, with a number of MSSM candidates in this fertile 
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patch of the landscape, it is now imperative to tackle the hard problems just mentioned. 
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A Physical states and string selection rules 

In this appendix, we discuss how to build consistent physical states. Furthermore, we list 
the string selection rules used in this work. Finally, we comment on an additional selection 
rule present in the literature: the 7 rule. We find that our construction of physical states 
is useful in order to apply the 7 rule correctly. It turns out that, in contrast to previous 
statements, the 7 rule does not further constrain allowed couplings. 

A.l Physical states 

An element of the space group g = {9^ ,naea) G S, where 9 is the twist and Bq, are the 
lattice basis vectors, corresponds to a boundary condition of a closed string [43,44]. For 
= (A; 7^ 0), the string is named untwisted string [twisted string). Focusing on its 
bosonic degrees of freedom in the six extra dimensions, the boundary condition reads 

A(r,a + 27r) = gX{T,a) , (A.l) 

where g is called the constructing element of the closed string. For each constructing 
element g, there exists a corresponding Hilbert space Tig of physical states. Using a 
mode expansion for A(r, cj), the general solutions of the string equation of motion with 
boundary condition Eq. (jA.ip can be written down. From these solutions one finds that 
twisted strings are localized at the fixed-point fg G IR^ corresponding to g (i.e. 9^fg + 
na^a = fg)- Furthermore, their quantization leads to the mass equation for left-movers 
(and the mass equation for right-movers is derived analogously). Focusing on the massless 
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case, the solutions are denoted bjo 

ksh>R ® \Pshh , (A.2) 

with shifted momenta qsh = Q+Vg and Pgh = P+Vg, where q and P he in the S0(8) weight 
lattice and Eg, x Eg, root lattice, respectively. The local twist and shift corresponding to 
the space group element g = {O^^naea) are defined by Vg = kv and Vg = kV + riaWa, 
respectively. Since the string is completely specified by its constructing element g and 
its left- and right-moving shifted momenta Pgh and Qs^, we write down a first ansatz for 
a physical state: 

Iphys) ~ ksh>R ® |^'sh)L ® b) (A.3) 

Up to now it is not guaranteed that a physical state is actually compatible with 
the orbifold. To ensure this compatibility, invariance of |phys) under the action of all 
elements of the orbifold group O C 5 G must be imposed (G is the embedding of 
S into the gauge degrees of freedom and is called the gauge twisting group). To do so, 
Eq. ()A.ip is multiplied by an arbitrary element h = {9 ^uiaea 

hX(T,a + 2-K) = hgX{T,a) (A.4) 
^ hX{T,a + 2tt) = hgh-^ hX{T,a) (A.5) 

Furthermore, the transformation properties of left- and right-movers under h are: 

ksh)R ® |^sh)L ksh)R ® \Psh)L , (A.6) 

wheri^ 

^ = g2vri [gsh-fh--Psh-Vh] (A 7) 

Now, we can distinguish two cases: 
Commuting elements: [h, g] = 

First, let us consider the transformation property of |phys) with respect to a commuting 
element h. In this case, Eq. (lATSj) yields 

hX{T,a + 2iT) = ghX{T,a), (A.8) 

i.e., the constructing element g is invariant under the action of h, 

\g) A \hgh-') = \g). (A.9) 



^^In this discussion we disregard oscillator states. Their inclusion is straightforward and does not 
change our conclusions. 

^^Here, we set "l>vac = 1 as discussed in [96]. 
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hX closes under the same constructing element g as X. Thus, both give rise to the same 
Hilbert space Tig 'H^g^-i = Tig- Furthermore, on the orbifold space M^/S" the string 
coordinates hX and X are identified. Thus, hX and X describe the same physical state. 

In summary, provided a constructing element g, we have shown that for commut- 
ing elements h, hX and X give rise to the same physical state from the same Hilbert 
space. Since h has to act as the identity on |phys), the following condition follows using 
Eqs. ([Oil . (IXej) and (TOj) : 

Qsh -Vh-Psh-Vh = mod 1 . (A.IO) 
Non— commuting elements: [h, g] ^ 

Next, considering a non-commuting element h in Eq. (|A.5p yields 

hX(T,a + 27r) = [hgh-^] hX{T,a) , (A.ll) 

i.e., the constructing element g is not invariant under the action of h, 

\g) ^ \hgh-^)^\g). (A.12) 

In the upstairs picture, i.e. in the covering space of the orbifold lR,^/5, one has 
different Hilbert spaces for the states with boundary conditions g and hgh^^. In this 
picture, Eq. (IXT2II says that h maps states from a given Hilbert space Tig onto a different 
Hilbert space Hhgh-i. Subsequent application of h then leads to the sequence 

^5 ^ '^hgh-^ ^ 'Hh2gh'2 Hh-igh-3 .... (A. 13) 

The crucial point is now that on the orbifold hX and X are identified. This means 
that, on the orbifold, the different Hilbert spaces TChngf^-n of the upstairs picture are 
to be combined into a single orbifold Hilbert space. Invariant states are then linear 
combinations of states from all H^n^^-n. Such linear combinations do, in general, involve 
relative phase factors (often called gamma-phase 7). So, the new ansatz for a physical 
state reads: 

Iphys) ~ ^(e-2"'"^ksh)R®|Psh>L®|/i"5^-")) 

n 

= \qsh)R0\Psh)L^(^e~^^'^^ \h^gh-^)^ , (A.14) 

where 7 = integer /A^, N being the order of the orbifold. The geometrical part of the 
linear combination transforms non-trivially under h 

^g-27rin7 ^ f^-n^ \ g2^i 7 ^ g-27ri n 7 . (A.15) 



^''Note that in all Tif^ngf^-n the left-moving momenta Psh of equivalent states are identical. The same 
holds for Qsh- 



29 



Since h has to act as the identity on |phys), the following condition follows using 
Eqs. (|A.6p . ()A.14|) and (jA.lSP for non-commuting elements: 

Qsh-Vh- Psh-Vh + 1 = modi. (A.16) 

Notice that 7 depends on h. Thus we can always choose j{h) such that this condition is 
satisfiecJ^. In principle, these steps have to be repeated for all non-commuting elements 
in order to ensure invariance of the physical state under the action of the whole orbifold 
group O C S G. The result for |phys) reads 

Iphys) = |(7sh)R®|^'sh)L® I Yl e-'^^'^^'^Uhgh-')] , (A.17) 

\h=l or [h,g]ytO ) 

where the summation over h is such that each term \hgh^^) appers only once. Note 
that the summation over h can be understood as a summation over all elements of the 
conjugacy class of g. 



Example 

To illustrate the construction of physical states, let us consider an example in the first 
twisted sector of the Zg-II orbifold. In the SU(3) lattice spanned by 63 and 64, there 
are three inequivalent fixed points associated to the constructing elements gi = {0, 0), 
92 = (^5 63) and g^ = {0, 63 + 64), or analogously gi = {0, ai 63 + hi 64) for i = 1, 2, 3 with 
ai = (0, 1, 1) and hi = (0, 0, 1). Then, restricting to the SU(3) lattice, the geometrical 
part of a physical state can be written as 

Y^^-2.i{n+mY< (n + m + ai) eg + (2m - n + hi) 64)) . (A.18) 

n,m 

Since the action of 9 in the SU(3) lattice has order 3, the only possible ^-eigenvalues of 
Eq. (jAlS]) have 7 = 0, ±\. In the case of 7 = 0, Eq. (jA.lSP is invariant under all rotations 
and translations for all three gi. However, if 7 = it^, the eigenvalue of Eq. (jA.lSP depends 
on gi: for the fixed point at the origin associated to gi, Eq. (lA.lSp is invariant under 0, 
but has an eigenvalue e^'^''^ ('^+') under (1, fees + ^64). Similarly, for the fixed points away 
from the origin, corresponding to gi {i 7^ 1), Eq. (jA.lSP picks up a phase q-'^'^^^ i'^i+bi) 
under (see Fig. [2|). It can be shown that for physical states 7 7^ is only possible in 
the presence of a Wilson line in the 63 and 64 directions. 



A. 2 String selection rules 

Consider the n-point correlation function of two fermions and n — 2 bosons [97, 98] 

(FFB...B). (A.19) 

""^^In this sense, building linear combinations and computing the 7 phase is not a projection condition. 
Note that 7(/i) is well-defined: if high^^ = h2gh2^ then 7(/ii) — 7(^2)- 
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Figure 2: Illustration of the 7-factor. The fixed point associated with the 
space group element 172 — 63) is invariant under (0, 63), but trans- 
forms into equivalent fixed points outside the fundamental domain under 
h = {9, 0). To form an eigenstate of {6, 0), one needs to build linear com- 
binations of the equivalent fixed points. The corresponding eigenvalues 
can be Ije^^^'''/^. 



The corresponding physical states shall be denoted by ^i, i = 1,... ,n. Then, in the 
field theory limit, a non- vanishing correlation function induces the following term in the 
superpotential 



A complete evaluation of Eq. (1A.19|1 has only been performed for 3-point couplings and 
yields a moduli dependent coupling strength [97,98,68,99]. 

On the other hand, symmetries of Eq. ()A.19P give rise to the so-called string se- 
lection rules. These rules determine whether a given coupling vanishes or not. We use 
the following notation: the constructing elements of are denoted hy gi £ S and their 
left- and right-moving shifted momenta, by Psh,i and q'sh,i) respectively. Then, the string 
selection rules read: 

1. Gauge invar iance 

The sum over all left-moving shifted momenta Psh,i must vanish: 



This translates to the field theoretic requirement of gauge invariance for allowed 
terms in the superpotential. 

2. Conservation of R charge 

R-charge is defined by 



n 



(A.20) 





Rt = Qs^,^-Nt + Nr fora = 0,...,3 
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where and A^*" are integer oscillator numbers, counting the number of excita- 
tions with oscillators d"" and q", respectively. Then the conditions [52] 

^Rf = modiV" for a = 1,2, 3 (A.23) 

i 

have to be imposed, where N'^ denotes the order of the twist component v'^ in 
the a-th complex plane, i.e. N°'v°' G Z (no summation). Here, two of the Ri come 
from fermions and the rest from bosons. For computational purposes, it is more 
convenient to use the purely bosonic notation, where Eq. ()A.23P becomes Yli RI = 
-1 mod N". 

This condition can be understood as a remnant of 10 dimensional Lorentz invari- 
ance. 

3. Space group selection rule 

The product of constructing elements gi must be the identity: 

n5. = (l,0). (A.24) 

i 

In terms of conjugate elements higih~^ of gi, this condition can be reformulated as 
W^Kgih-^ = {l,v) with V G - ^'^'OA [100]. 

This selection rule can be visualized as the geometrical ability of twisted strings to 
join. 



A. 3 On the need for a 7 selection rule 

In the literature, there exists an additional selection rule, here referred to as the 7 rule. 
In our notation, it reads [52,68] 

^ 7i = mod 1 , (A.25) 

i 

where 7^ denotes the gamma-phase of ^'j. In this section, we argue that, in contrast to 
previous statements, a fully consistent approach yields to automatic fulfillment of the 7 
rule. 

The correlation function corresponding to the coupling 
^1^2...^n (A.26) 

should be invariant under the action of the full space group. Let us assume first that 
the states corresponded to linear combinations of equivalent fixed points within the 
fundamental domain of the torus (see e.g. [52,68, 101]). For example, in the case of the 
Zg-II orbifold only fixed points in the G2 lattice could form linear combinations. Under 
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this assumption, different states would be eigenstates witli respect to different space 
group elements. So one could not transform the coupling Eq. (|A.26j) with a given h = 
{6\maea)- Thus the fully consistent approach for building invariant linear combinations, 
as presented in Appendix lA.ll is necessary. In this case, we can compute the gamma- 
phase for all states ^'j from Eq. ()A.16p . i.e. 7^ = 7i(/i) for arbitrary h = {6\maea)- But 
since allowed couplings already fulfill the selection rules Eqs. ()A.2ip and (|A.23p . the 7 
rule is satisfied triviallj0: 

7i(/i) = Psh,i ■ Vh - qsh,i • Vh , (A.27) 




(A.28) 



= mod 1 . 

Thus, the 7 rule in the fully consistent approach is not a selection rule. It is a consequence 
of other selection rules and invariance of the states. We therefore conclude that the 
coupling must only satisfy gauge invariance, i?-charge conservation and the space group 
selection rule. 

This has important consequences. For example, in the model Al of [52], there is 
no mass term for the exotics 52^2 up to order 9 in singlets. However, we find that the 
coupling q2q2 39815822333 is allowed by the selection rules of Appendix lA.2[ Further, using 
the prescription of Appendix lA.H the gamma-phases of the corresponding physical states 
are ji = (|, 0, 0, |, |, 0) for h = {9, 0) , which sum up to 2. This is in contrast to [52], 
where 7^ = (0, 0, 0, ^, |, 0) and linear combinations were built differently. 



B D— flatness 



In this appendix, a simple method is reviewed that allows to analyze L>-flatness. It also 
provides a simple test whether it is possible to cancel the FI term with a given set of 
fields. 

Let us start by briefly reviewing the issue of D-flatness and cancellation of the FI 
term [102-106]. In supersymmetric theories, there is the so-called D-term potential. In 
the case of a U(l) gauge theory it is given by 

2 



Vd oc 



(B.l) 



Consider as a first example a U(l) gauge theory with two fields (f)± carrying the 
charges ±1. Clearly, as long as \4>+\ = \4>-\, Vd vanishes. That is, one has a D-Rat 
direction, parametrized by x = |(/)+| = |(/>_|. 



^Also in the presence of oscillators, the 7 rule is satisfied automatically. 
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Consider now a theory with one field with charge 2 and two fields {4>2,4>3) 
with charges —10 Then we have many flat directions, described by the roots of the 
equation 2|(/)ip — \(j)2\'^ — I'Psl'^ = 0. It is convenient to associate these directions to the 
(holomorphic) monomials 



h (Pl , </>l 03 > 01 02 03 



respectively. That is, a monomial 0"^ 02^ • • • 0^*" represents a flat direction, defined by 
the relation 



I0A 



and |0j| = for Uj = . 



The crucial feature of such monomials is that they are (obviously) gauge invariant. More 
precisely, every holomorphic gauge invariant monomial represents a D-flat direction [102]. 

It is, however, clear that there is only a finite number of linearly independent D-flat 
directions. In the previous example, the third direction is not independent of the other 
two. In other words, the requirement Vd = poses only one constraint on the three real 
variables (|0jp) entering (|B.ip . The space of absolute values |0i| is 2-dimensional. The 
power of using the monomials is that checking whether certain monomials are linearly 
independent or not is fairly simple: identify with each monomial the vector of exponents, 
V = (ni,n2 . . . ). The directions are independent if and only if the vectors are linearly 
independent. In the previous example one would get the vectors (1,2,0), (1,0,2), and 
(1, 1, 1), out of which only two are linearly independent. 

It is also clear how to obtain these vectors: all of them are orthogonal to the vector of 
charges q = {qi,q2 ■ ■ ■)■ That is, the problem of finding the above monomials (and thus 
the D-Rat directions) is reduced to the problem of finding vectors v with the following 
properties: 

1. q-v = 0, 

2. Vi £ No. 

The property that the Vi be integer-valued does not pose a constraint in our models: 
since the charges are rational, one can rescale any v having the first property such as to 
have integer entries. However, the requirement that the entries be non-negative, which 
reflects that the monomials ought to be holomorphic, is a constraint. 

The discussion so far can easily be extended to U(l)'^ theories. Here the Z)-term 
potential is 

2 



E 



,0) 



(B.2) 



^''For the moment, we ignore anomalies. 
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where qf^ is the charge of the field (f)i under the j*^ U(l) factor. Now a D-flat direction has 
to satisfy the above constraints for each U(l) factor separately. Again, it is advantageous 
to represent D-flat directions by holomorphic gauge invariant monomials (dubbed 'HIMs' 
in the literature [106]). Then the vector v of exponents has to be orthogonal to every 



charge vector q^^^ 
matrix 



). In other words, v has to be in the kernel of the charge 



( <A 



(1) 



Q ■ V = , with Q 



(2) 
9l 



(n) 



4" 



(n) 



(B.3) 



... J 



Hence, the problem of finding the D-fiat directions of a U(l)" gauge theory is reduced 
to the task of calculating the kernel of the charge matrix Q, and to forming linear 
combinations of elements of this kernel in such a way that the entries are non-negative 
integers. The maximal linear independent set of such linear combinations is in one-to-one 
correspondence with the independent D-fiat directions. 

Next, let us comment on what happens if there are non-Abelian gauge factors. Then 
the D-term potential is to be amended by 

2 



T/non— Abclian 



E 



(B.4) 



with Tq denoting the group generators. It is straightforward to see that the results 
obtained so far generalize to the non-Abelian case [102]: the D-Qat directions are again 
in correspondence with holomorphic gauge invariant monomials. That is, one can amend 
the monomials discussed so far such as to include fields transforming non-trivially under 
non-Abelian gauge factors, as long as these fields are contracted in such a way that the 
monomials are gauge invariant. 

Finally, let us review the issue of cancelling the FI term. For an 'anomalous' U(l), 
the D-term potential ()B.ip gets modified to 

2 



Vi 



D 



oc 



(B.5) 



where in our convention ^ > 0. To cancel the FI term one thus has to find a holomorphic 
monomial, 



with net negative charge under U(l) 



1, i-e. 



Uiq™ < 0. 



(B.6) 



(B.7) 
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To summarize, the D-Rat directions are in one-to-one correspondence with holomor- 
phic gauge invariant monomials. In the Abehan case, such monomials can be identified 
with elements of the kernel of the charge matrix Q with non-negative integer entries. 
Cancelation of the FI term requires the existence of a holomorphic monomial with net 
negative charge under U(l)anoirn which is gauge invariant with respect to all other group 
factors. 



C Family reflection symmetry and Matter Parity, 

We would like to define an effective low energy theory which preserves R parity. This has 
the advantage of greatly reducing the number of arbitrary parameters in the superpo- 
tential, forbidding dimension 3 and 4 baryon or lepton number violating operators, and 
preserving a viable dark matter candidate, i.e. the LSP. Our strategy for accomplishing 
this is, in principle, quite simple. We make use of "family reflection symmetry" or "mat- 
ter parity" defined as a discrete subgroup of U(1)b-l- This is a global symmetry 
(commuting with supersymmetry) which is even on the Higgs doublets and odd on all 
SM quark and lepton fields. It forbids the following dangerous baryon or lepton number 
violating operators, 

udd, qde , e£e and ihu- (C.l) 

On the other hand, it allows quark and lepton Yukawa couplings as well as the Major ana 
neutrino mass operator Pz^. 

Consider the effective operators 

0{si... Sn) with O in jUI]) , (C.2) 
VV{s',...s'J . (C.3) 

We want to forbid the dangerous proton decay operators (lC.2p . while allowing for Ma- 
jor ana neutrino masses ()C.3|) . This puts a constraint on the B — L charges of the SM 
singlets which get non-zero VEVs. In particular, it requires 

-1 + J^g, / (C.4) 

i 

for any set of singlets with non-zero VEVs, where qi are the B ~ L charges. In addition, 

i 

must be satisfied for at least one singlet configuration with B — L charges q'j^. In our 
theory, the singlets come in pairs with opposite B — L charges and these charges are 
rational. Then the relevant solution to the above equations is 

^ 2li + l ^ ^ 
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quarks 


leptons 


K 


hd 


udd 


qM 


Ue 




B-L 


1/3 


-1 








-1 


-1 


-1 


-1 


R 


gia7r/3 




1 


1 











Table A: B — L and i? charges for SM particles, with opposite B — L 
charges for anti-particles, and for baryon and lepton number violating 
operators with dimension < 4. 



for integer ki and h, with the additional condition that adding/subtracting the numera- 
tors of Qi can yield 2, i.e. 

^hNi = 1, (C.7) 

i 

for some Ni € Z. For example, the numerators can differ by 2. If there are only two fields 
with charges ±qa, the corresponding constraint is qa = ±2/(2/ + 1). 

The above singlet VEVs break U(1)b_l to a discrete subgroup. Consider an element 
of U(1)b_l defined by 

R{a) = e'^"'^-^ . (C.8) 

The B — L and R{a) charges of SM particles are given in Table [H The choice a = 3 
corresponds to family reflection symmetry (FRS). If only the singlets satisfying 

qs-d^ = (C.9) 

obtain VEVs and there is at least one singlet for which Z 7^ 0, then \J{1) b-l is broken 
to R{3) = FRS. Clearly, products of these singlets can contribute to Yukawa couplings 
for quarks and leptons. Further, products of singlets with B — L charge —2/3 (or, more 
generally, those with B — L charge 2/3 and —4/3, etc.) can generate Majorana neutrino 
masses. On the other hand, the proton decay operators are forbidden. 

It is possible to generalize FRS to a Zat group. In general as long as a 7^ 2Z, R{a) 
will forbid the dangerous operators, Eq. (jC.ip . and allow all Yukawa couplings. Consider 
a field (j) C {s] with qs-i^'P) = / • Such a field breaks V{1)b-l to the subgroup 
R{a) with a = 2//. The effective Majorana neutrino mass operator vUcj)"' is allowed for 
[2 + nf]/ f = Z or a = 2// = Z. Hence, for odd a, we can both forbid the dangerous 
operators, Eq. (jC.ip and obtain non-zero Majorana neutrino mass. The corresponding 
constraint on / is then / = ib2/(2Z + 1), as expected. For example, 

1. / = ±2, a = 1 gives R = e'^^^^-L g 

2. / = ±2/3, a = 3 gives R = e^i'^es-L ^ 
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3. / = ±2/5, a = 5 gives R = e^'^^^^-L ^ ^q, 

4. / = ±2/7, a = 7 gives R = c'^'^'Ib-l ^ Zq, 

5. / = ±2/9, a = 9 gives R = e^'^^^^-L ^ ^2. 

This is easily generalized to configurations with many different singlets getting VEVs 
(with the constraints given in (IC.6P and (IC.7P ). These conserve matter parity, ^2^. 



D Search for B — L and R parity 

Our search for \J{1) b~l is based on the methods developed in ref. [62] (for an earlier 
discussion of \J{1) b-l and its applications see [53] and [40]). In Tab. [H we list the 
standard model particle content with their hypercharge and B — L charges. 



q 


(3,2) 


1/6, 1/3 




(1,2) 


-1/2,-1 


K 


(1,2) 


1/2, 


u 


(3,1) 


-2/3,-1/3 


e 


(1,1) 


1, 1 




(1,2) 


-1/2, 


d 


(3,1) 


1/3,-1/3 


U 


(1,1) 


0, 1 





Table 5: Matter content of the standard model, where the sub- 
scripts denote hypercharge and B — L, respectively. In our conventions, 

Q = T3L + Y. 

The choice for \]{1)b-l depends on the choice of hypercharge in the first place. In 
this publication, we do not take the most general approach, but assume that hypercharge 
is given by SO(IO) D SU(5) D SU(3) x SU(2) x U(l)y. Furthermore, we demand that the 
first and second families come from 16-plets localized in the first twisted sector, whereas 
the multiplets of the third family may come from any sector of the theory. 

To find a suitable \J{1) b-l, we proceed as follows. In general, the shift and Wilson 
lines break the gauge group in 10 dimensions 

Eg X non-Abelian x U(l)" . (D.l) 

The U(l) generators are n linearly independent directions U in the root lattice of Eg x Eg 
that are orthogonal to the simple roots of the unbroken non-abelian gauge group. For 
the B — L direction, we make the general ansatz 

Ib-l = xi ti ± X2t2 ± . . . ± a;„t„ . (D.2) 

The B—L charge of a particular representation is given by the scalar product of its highest 
weight and B — L. We denote the highest weights of the left-handed quark doublets and 
of the right-handed quark singlets by Aj, i = 1, . . . , 9. Note that the first two families are 
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fixed, and we loop over all representations which have the right quantum numbers to be 
the quarks of the third generation. For each such choice, we have 



Aj • tB-L = Xi Aj • ti + X2 Aj • t2 + . . . + x„ Aj • tn , (D.3) 

which is a system of 9 linear equations. Although one may think that these 9 equations 
severely constrain the values of Xj, this is not true. In general, the system will be under- 
determined, since the quarks may differ by localization, but not necessarily by the highest 
weights of their gauge representations. In order to account for the B — L charges of the 
leptons and Higgses, and for the absence of chiral exotics, we set up necessary, but in 
general not sufficient, linear constraints: 

^ qB-L = 1 , ^ qs-L = - 2 , 

(3,2), (3,2) (3,1), (3,1) 

Y^qB.L= -3, 1B-L= 3. (D.4) 

(1,2) (1,1), y^o 

For readability, we use qB-L in the above equations as a shorthand for Afc-t^-L, where A^ 
runs over the highest weights of the representations in the sum. Note that the sum over 
e.g. the (3, 1) and (3, 1) representations reduces to that over the right-handed quarks u 
and d alone, since we assume pairs of exotic particles to carry B — L charge assignments 
that are equal in magnitude but opposite in sign. Another linear constraint comes from 
the requirement that U(1)b-l be non-anomalous. If the model has an anomalous U(l) 
direction that we will denote by tanomi we demand that it be orthogonal to B — L: 

(Xi ti + X2 t2 + . . . + Xn t„) • tanom = . (D.5) 

In principle, we could write down many more linear conditions, e.g. we could demand 
that the sum of the B — L charges of u-type quarks like (3, l)_2/3,-i/3 and (3, 1)2/3,1/3 
alone gives -1, which is more constraining than Eq. ()D.4p . However, our experience shows 
that it is more practical to drop these conditions, since the non-linear equations to be 
introduced below are constraining enough. 

The linear equations by no means exhaust the constraints we may require to be 
fulfilled by a vector-like spectrum. In particular, there are the cubic, 

(3,2), (3,2) (3,1), (3,1) 

Y^{qB-Lf=-S, Yl iQB-Lf=3, (D.6) 

(1,2) (1,1), Y^O 

and the quintic, 

E_ ('?^-^)' = ^ ' E_ = - 1[ ' 

(3,2), (3,2) (3,1), (3,1) 
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Y^{qs.Lf= -3, (1B-Lf= 3, (D.7) 

(1,2) (1,1), y^o 

constraints. As before, the sum over all (1,1) representations with non- vanishing hyper- 
charge reduces to that of the right-handed electrons that carry B — L charge +1, since the 
exotic particles come in vector-like pairs so that their contribution to the sum vanishes. 

This leaves us with a set of highly non-linear equations to be solved. To this end, 
we used the computer algebra system Singular [107]. The following cases need to be 
distinguished, (i) The number of solutions is finite, (ii) The solutions are given by con- 
tinuous parameters, and the relations intertwining these parameters are linear, (iii) The 
solutions are given by continuous parameters, but this time, the relations intertwining 
the parameters are non-linear. B — L directions which lead to irrational charges for the 
exotics are discarded in all three cases. 

In the first case, we calculate the spectra for the B — L generators and check that they 
arc really vcctor-likc. A given model is particularly interesting if its spectrum contains 
standard model singlets with B — L charges 2n/ (2m + 1) for n, m € N whose VEVs break 
the continuous \]{1)b-l symmetry to a discrete subgroup which may play the role of a 
generalized R parity. In the second case, we use the parametric freedom we have at hand 
to assign B — L charges of 2n/ (2m -|- 1) for m, n = 0, 1, 2, 3 in all possible combinations 
to standard model singlets and thus generate a list of B — L directions. In the third case, 
we specialize to a numerical value and check whether the spectrum satisfies our criteria. 
We do not perform a complete search, because these cases are rare and the analysis more 
complicated, thus not very rewarding. 



E Details of Model 1 

The model is defined by the shifts and Wilson lines 

/111 \ /I 1 1 1 1 1 1 1\ , , 

^ = (3'-2'-2'«'°'«'°'°j U'-6'-2'-2'-2'-2'-2'2j ' ^^"^^^ 

/llll \/ 7 79 7\ 

= (o,-5,--.--,5,0,0,oj (^4, -3, --,-4, -3, --,--.-] , (E.lb) 

A possible second order 2 Wilson line is set to zero. The gauge group after compactifi- 
cation is 

Gorbifoid = SU(3) X SU(2) X [SU(4) x SU(2)'] x U(l)9 . (E.2) 
The U(l) generators can be chosen as 
ti = ty = (0,0,0,i,i,-|,-i,-|) (0,0,0,0,0,0,0,0) , (E.Sa) 
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t2 = 


(1,0,0,0,0,0,0,0) 


(0,0,0,0,0,0,0,0) , 


(E.3b) 


t3 = 


(0,1,0,0,0,0,0,0) 


(0,0,0,0,0,0,0,0) , 


(E.3c) 


U = 


(0,0,1,0,0,0,0,0) 


(0,0,0,0,0,0,0,0) , 


(E.3d) 


t5 = 


(0,0,0,1,1,1,1,1) 


(0,0,0,0,0,0,0,0) , 


(E.3e) 


t6 = 


(0,0,0,0,0,0,0,0) 


(0,1,0,0,0,0,0,0) , 


(E.3f) 


t7 = 


(0,0,0,0,0,0,0,0) 


(-1,0,0,1,0,0,0,0) , 


(E.3g) 


ts = 


(0,0,0,0,0,0,0,0) 


(0,0,0,0,1,0,0,0) , 


(E.3h) 


tg = 


(0,0,0,0,0,0,0,0) 


(0,0,0,0,0,1,0,0) . 


(E.3i) 



The 'anomalous' U(l) is generated by 

9 



^Oitj , where {aj 



i=l 



/n 2 Q _5 1 _i 1 2 i\ 

<(^u, 3, u, 3,3, 3, 3, ^, 3 j> 



(E.4) 



The sum of anomalous charges is 
296 

tr tanom — ^ . 



(E.5) 



E.l Spectrum 



Table 6: The spectrum of model 1 in terms of left-chiral states. 
The U(l) charges refer to the basis of generators (IE.3j) . Ri denote 
i?-charges. Note, the symbol * means the state is not localized in 
the corresponding directions. 
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E.2 Mass matrices for exotics in Model lA 

Here we show that that all exotics can be made massive. The exotic's mass terms are 

Xi {Mxx)ij Xj . (E.6) 
In the following, we list the structure of the corresponding mass matrices. 
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E.3 Mass matrices for exotics in Model IB 

Here we show that that all exotics can be made massive. The exotic's mass terms are 
Xi {Mxx)ij Xj . (E.8) 
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In the following, we list the structure of the corresponding mass matrices. 
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(E.9h) 
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F Details of Model 2 

The model is defined by the shift and Wilson lines [40] 



W2 = 



1 


1 


1 


3' 


2' 


2' 


1 


1 


1 


4' 


4' 


4' 



1 1 1 1 1 1 1~ 

2' 2'6'6'6'6'6'6, 



1 

2 


1 

2 


1 
2 


1 

2 


2' 2/' 


1, 


5 


3 


1 


5 3 


2' 


2' 


2' 


2' 2' 


6, 


7 
3' 


4 
3' 


-5, 


-3, 3) . 



1 1 1 IWi _-| _5 _3 _1 _5 _3 3\ 
4) 4'4' 4/ V ' ' 2' 2' 2' 2' 2' 2/ 



(F.la) 
(F.lb) 
(F.lc) 



We use the U(l) generators 



tl 


= ty = (0,0,0,i,i,-i,-i,-i) (0,0,0,0,0,0,0,0) 


(F.2a) 


t2 


= (1, 0, 0, 0, 0, 0, 0, 0) (0, 0, 0, 0, 0, 0, 0, 0) , 


(F.2b) 


t3 


= (0, 1, 0, 0, 0, 0, 0, 0) (0, 0, 0, 0, 0, 0, 0, 0) , 


(F.2c) 


t4 


= (0, 0, 1, 0, 0, 0, 0, 0) (0, 0, 0, 0, 0, 0, 0, 0) , 


(F.2d) 


t5 


= (0,0, 0,1, 1,1, 1,1) (0,0, 0,0, 0,0, 0,0) , 


(F.2e) 


t6 


= (0, 0, 0, 0, 0, 0, 0, 0) (1, 0, 0, 0, 0, 0, 0, 0) , 


(F.2f) 


t7 


= (0, 0, 0, 0, 0, 0, 0, 0) (0, 1, 0, 0, 0, 0, 0, 0) , 


(F.2g) 


t8 


= (0, 0, 0, 0, 0, 0, 0, 0) (0, 0, 0, 1, 1, 0, 0, 0) . 


(F.2h) 



The anomalous U(l) generator can be expressed through a linear combination of all U(l) 
generators, 

5 1 2 2 2' 



^Qtj, where q = ^0,^,-1, 



3'3'3' 3' 3 



The sum of anomalous charges is 

416 
~3~ 



trta 
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(F.3) 



(F.4) 



F.l Spectrum 
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Table 7: The spectrum of model 2 in terms of left-chiral states. 
The U(l) charges refer to the basis of generators (|F.2p . Ri denote 
i?-charges. Note, the symbol * means the state is not localized in 
the corresponding directions. 
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